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Abstract 

We construct an approximate scheme based on the concept of the spontaneous sym- 
metry breakdown, satisfying the Goldstone theorem, for finite volume Bose-Einstein 
condensed gases in both zero and finite temperature cases. In this paper, we discuss 
the Bose-Einstein condensation in a box with periodic boundary condition and do 
not assume the thermodynamic limit. When energy spectrum is discrete, we found 
that it is necessary to deal with the Nambu-Goldstone mode explicitly without the 
Bogoliubov's prescription, in which zero-mode creation- and annihilation-operators 
are replaced with a c-number by hand, for satisfying the Goldstone theorem. Fur- 
thermore, we confirm that the unitarily inequivalence of vacua in the spontaneous 
symmetry breakdown is true for the finite volume system. 
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1 Introduction 



Various theoretical studies on weakly interacting dilute Bose-Einstein con- 
densed systems have been done over many years. The excitation spectrum in a 
homogeneous system was first obtained by Bogoliubov in 1947 [1]. Afterwards 
the quantum correction in zero-temperature case was evaluated by Hugenholtz 
and Pines [2] . In course of their discussion, they showed the Hugenholtz-Pines 
(HP) theorem, which relates the chemical potential with the self-energy in a 
nonperturbative way. In these works one follows the Bogoliubov's prescription 
[1], in which the zero-mode creation- and annihilation-operators do and CLq £1X6 
replaced with a c-number by hand, where N c is the number of conden- 
sate particles, and the c-number is regarded as the order parameter. Though 
this prescription was simple and successful in deriving the excitation spec- 
trum, its naive use is not consistent in quantum theory because the canonical 
commutation relations (CCRs) are not respected [3]. 

In modern quantum field theory, the order parameter is introduced as a vac- 
uum expectation value of a quantum field. The self-consistent mechanism for 
creating the order parameter is the spontaneous symmetry breakdown (SSB). 
We have many examples of the SSB, appearing in ferromagnetic, supercon- 
ducting, crystal orders, and so on [4]. The scenario of the SSB reads as follows: 
Suppose that an action or a Lagrangian is invariant under a continuous trans- 
formation but the vacuum does not show the original symmetry. Then the 
non-zero vacuum expectation value of the quantum field is induced, and is 
identified as the order parameter. In general, when a continuous symmetry 
is broken spontaneously, there must exist a gapless mode, called the Nambu- 
Goldstone (NG) mode, which is necessary to keep the original invariance of 
the action. This statement is the Goldstone theorem [5] , and it is proven from 
the Ward-Takahashi (WT) relations [6,4]. The WT relations are in turn iden- 
tities, derived from the CCRs, the Heisenberg equation of motion, and the 
transformation property. 

The dilute gas systems of bosonic neutral atoms at recent experiments [7,8,9], 
forming Bose-Einstein condensates (BECs), are described by the model La- 
grangian invariant under a global phase transformation. The appearance of 
BECs can be interpreted as a spontaneous breakdown of the global phase 
symmetry. For such a model in homogeneous situations, Hohenberg and Mar- 
tin [10] showed that the HP theorem is a consequence of the WT relation 
with respect to the global phase transformation first, and generalized it to a 
finite temperature case. Recently, Boyanovsky et al [11] showed that the HP 
theorem holds at one-loop level in a homogeneous BEC system. 

The concept of the SSB is inherent in quantum field theory but not in quantum 
mechanics because it is closely related to an infinite number of degrees of free- 
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dom, causing many unitarily inequivalent vacua. For spatially homogeneous 
systems of quantum field one finds inequivalent vacua under the thermody- 
namics limit, defined as that of infinite volume while its ratio to the number 
of particles is kept finite [4]. 

The whole story above applies in the case of homogeneous systems and can 
not be extended to spatially inhomogeneous cases without due consideration. 
The BEC systems of neutral atoms at recent experiments have finite system 
size and spatial inhomogeneity due to trapping potentials. For such systems 
of finite size, the invariance of spatial translation is lost and the momentum 
is no longer a good quantum number, and the energy spectrum including 
the NG mode becomes discrete. The issue of existence of inequivalent vacua 
for such finite volume systems has not been considered thoroughly, so the 
occurrence of the SSB in the systems is not trivial. Meanwhile, it has been ob- 
served in the recent experiment on trapped Bose gases [12] that a condensate 
has a fixed phase. Furthermore, it is claimed that the Bogoliubov spectrum 
[13] and the Bogoliubov transformation [14] were observed in the two-photon 
Bragg scattering process. Even in the case of discrete spectrum, the theoreti- 
cal formulation of quantum field with introduction of the order parameter is 
consistent with the experiment [15]. These experiments seem to support that 
as a result of the SSB the order parameter exists with a fixed phase, not only 
in homogeneous infinite volume system, but also in inhomogeneous system of 
finite size. 

When an order parameter appears, a single vacuum should be selected among 
infinite possible ones and a quasi-particle picture should be established in a 
manner consistent with the Goldstone theorem. This is a vital but not trivial 
step in theoretical calculation. One needs to invent a selection procedure in 
theory. In the Bogoliubov prescription for a homogeneous system, the vacuum 
is selected through replacing a and aj, with V^c by hand, and the Bogoliubov 
transformation of p ^ modes diagonalizes the unperturbed Hamiltonian. It 
is remarked that the zero-mode operators are absent in the formulation so 
that the CCRs for the field operators and the Goldstone theorem can not 
be satisfied exactly. But since the point of the zero-energy is embedded in a 
continuous spectrum, this fact is easily overlooked. In case of finite size sys- 
tem such as the trapped BEC systems with discrete spectrum, the Bogoliubov 
prescription clearly fails: The CCRs of the quantum field and the Goldstone 
theorem are violated due to the absence of the discrete zero-mode. Conven- 
tionally, in the SSB of quantum field theory, one adopts the systematic method 
to introduce an infinitesimal symmetry breaking term, known as Bogoliubov's 
quasi- average [16], and one vacuum is selected among many. Note that the 
zero-mode is naturally included in this method in the symmetric limit. We 
have proposed the way of introducing the infinitesimal symmetry breaking 
term for the trapped BEC system and have explicitly shown that the WT 
relations and the Goldstone theorem hold at tree level [3]. It is also shown 
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that the vacua belonging to the order parameters with different phases are 
unitarily inequivalent to each other even for the trapped BEC [17]. 

Our goal is a formulation of a systematic and consistent treatment of inhomo- 
geneous quantum field system of finite size, corresponding to the recent ex- 
periments of the trapped BEC. We mean by the word "consistent treatment" 
that the CCR, the WT relations and the Goldstone theorem are respected. 
As mentioned above, we have studied this subject at tree level [3], and should 
extend our discussions to any loop level. At present, such extension is not easy 
due to the complex structure of the unperturbed propagators. 

Instead, in this paper, we consider a BEC system, without trapping poten- 
tial but confined in a box with periodic boundary condition. We keep a finite 
volume size V. To select a vacuum, we introduce an infinitesimal symmetry 
breaking term, characterized by an infinitesimal parameter e. What is impor- 
tant in our present study is that the limit e — > is taken but the limit V — > oo 
is not, namely that the thermodynamic limit is not considered. We will see 
that inequivalent vacua emerge in the limit e — > with finite V (without 
the thermodynamic limit). While the energy spectrum is discrete, the Fourier 
(momentum) representation is available in this model, which makes calcula- 
tions of loop expansions tractable. It is expected that a model with trapping 
potential mostly shares essential theoretical features with this model. 

Explicitly we show that the HP theorem is derived from the WT relations even 
in this finite volume system. The quasi-particle picture in which the zero-mode 
(NG mode) is present is constructed at tree and one-loop levels. To extend 
the results at zero temperature to those at finite one, we employ thermofield 
Dynamics (TFD) [4,18] which is a real time canonical formalism of thermal 
field theory with doubled degrees of freedom. 

This paper is organized as follows. In Section. 2, we give a formulation of the 
model Lagrangian density with an infinitesimal symmetry breaking term at 
zero temperature, and construct the quasi-particle picture consistently with 
the CCRs and the Goldstone theorem. The order parameter is introduced 
without using the Bogoliubov's prescription. We give the quantum correction 
which keeps the HP theorem without taking the infinite volume limit. We also 
show that the vacuum in a broken phase is orthogonal to that in as symmetric 
phase. In Section. 3, we use TFD to extend the results in Section. 2 to finite 
temperature case. We review the TFD formalism briefly and investigate the 
WT relations at each loop level similarly to the zero temperature case. In this 
case, we also derive the HP theorem from the WT relations at arbitrary loop 
level. Section 4 is devoted to a summary and conclusion. In Appendix, we 
derive the relation between the vcua in the broken and symmetric phase. 
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2 Zero Temperature Case 



We formulate a field-theoretical treatment for the BECs in a box with periodic 
boundary condition on the field variables. In this section, we discuss the zero 
temperature case. First we confine our discussion to tree level, and will see 
that the WT relation is satisfied and that the unitarily inequivalent vacua 
are realized. Next we derive the HP theorem from the WT relations at any 
loop level and show explicitly that our one-loop calculation satisfies the HP 
theorem. 



2. 1 Model Lagrangian density and Hamiltonian 



We start with the following Lagrangian density, which describes a weakly 
interacting Bose gas in a box whose volume is denoted by V, 



(1) 



Here, to, \i and g are the mass of a boson, the chemical potential and the 
coupling constant, respectively, and x stands for the space-time coordinate 
(x,t). The periodic boundary condition in each direction is imposed on the 
field variables: 

V(x + L,y,z,t) = V(x,y,z,t) 

V(x,y + L,z,t) = V(x,y,z,t) , (2) 
V(x,y,z + L,t) = V(x,y,z,t) 



where L being the length of one side of the cube, V = L 3 . It is easy to see that 
this Lagrangian density is invariant under the global phase transformation: 



¥(x)^e- w ^(x), (3) 

where 6 is a real constant. When a uniform BEC is created, the global phase 
symmetry is spontaneously broken, and then the quantum field is divided 
into a classical constant field v and a quantum field (p(x) in the terminology 
of the canonical operator formalism: 

^(x) = v + (p(x) . (4) 
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We take a real value for v, which does not affect generality of the following 
discussions. The classical field v is called the order parameter and is also 
expressed as v 2 = n c in terms of the density of condensate particles n c . The 
classical field v may be defined as an expectation value of the Heisenberg field 
with respect to the vacuum \Q), 

(n\v(x)\n) = v, (5) 

or equivalently 

(n\ip(x)\n) = o. (6) 

Let us introduce an artificial symmetry breaking term 

C £ = (se)v \v(x) + &(x)] , (7) 



and add it to the original Lagrangian density (1), then the total Lagrangian 
density C tot becomes 

C tot = C + C £ . (8) 



Here, e is an infinitesimal dimensionless parameter and e represents a typical 
energy scale of the system. Obviously the total Lagrangian C tot is not invariant 
under the global phase transformation (3). The parameter e is taken to be 
vanishing at the final stage of the calculation, so that the original symmetry 
is restored. 

We move to the canonical formalism in the interaction representation. The 
canonical conjugate of <p(x) is ir(x) = iTup^x). Then the CCRs are as follows: 

[ V (x,t),^(x',t)] = 5 3 (x-x'), (9) 

[<p(x, t), <p{x', t)\ = [<p\x, t), <p\x', t)] = 0. (10) 

Now, we have the total Hamiltonian 

H = H + H int + const. , (11) 

where the unperturbed Hamiltonian H is given as 



Ho = Jd 3 x ^( x )(-^-V 2 -f?j<p(x) 
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\-^— ^A(p\x)(p(x) + (p(x)tp(x) + <p*{x)<p\x)} 



(12) 



and the perturbative Hamiltonian H int is defined as 



H- mk = v(-n + gv 2 - ee) J d 3 x {<p(x) + v\ x )} 

+ gv J d 3 x (p J< (x)(p J< (x)(p(x) + (p^ (x)(p(x)(f(x) 



(13) 



It is convenient to define the 2x2- matrix propagator for the quantum fields 
(p(x) and its Hermite conjugate (p^(x) by 



G(x - x') 



Gn(x-x') G 12 (x-x') 
G 21 (x-x') G 22 (x-x') 

-i(n\T[ip(x)^(x')}\n) -i(n\T[ip(x)ip(x')]\n) 

-i(Q\T[^{x)^(x')]\Q) -i(Q\T[^{x) V {x')]\Q) 



(14) 



where T is the symbol for time-ordered product. Its Fourier transform is given 
by 



G( P ) 



Gn( P ) G 12 (p) 
G 21 (p) G 22 {p) 



d 4 x 
(2nTi) [ 



■G(x)e^ p - x -° Jt) , 



(15) 



with the notation of p — (ou,p). It is mentioned that one has the following 
relations among the matrix elements: 



G n (p) = G 22 (-p) , G 12 (p) = G 21 {p) . 



(16) 



2.2 Derivation of the Ward-Takahashi relation 



It is well-known that the Goldstone theorem follows from the WT relations. 
We adopt the loop expansion, because the WT relations hold at each loop 
level in general. 
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First we review the WT relation in general case. Consider an infinitesimal 
transformation: 

-> V'(x) = + (f 7) 



with the infinitesimal parameter £. The change in the Lagrangian density is 
denoted by 

£5£ = C[V(x)] - £[*(*)] . (18) 



The Nother theorem implies 

6£ = & i N ll , (19) 



where 



N " = • (20) 



Now, we define the Nother charge 

N(t) = ^Jd 3 xN (x), (21) 

which generates the transformation in the commutation relation 

[*(x),N(t)] tm = t = m(x). (22) 

Then, the Nother theorem (19) leads to 

N(t) = ^Jd 3 x5£. (23) 

One can easily derive the following relation, 

^(n\T[N(t)*( Xl )...*(x n )]\n) 

= E - t a ){n\T[*( Xl ) ■ ■ ■ [N(t), *(x a )} ■ ■ ■ *(x n )]\Sl) 

o=l 

+ (n\T[N(t)*( Xl )..-*(x n )]\n). (24) 
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Integrating both sides over the entire time domain and dropping the surface 
term at t — ±00, we obtain 

n „ 

£ ifc(ft|T[tt(a;i) • --6*(x a ) ■ ■■^(x n )]\Q) = / d 3 x(n\T[5C(x)^( Xl ) ■ ■ • *(x n )]|fi) .(25) 
0=1 ^ 



This relation, derived independently of any approximate scheme, is called the 
WT relation. 

In this paper, we are interested in the following infinitesimal global phase 
transformation: 

8V(x) = W(x) . (26) 



The Nother charge is given by 

N(t) = - Jd 3 x¥(x)^(x) , (27) 



indeed [N (t) , ^ (x)] tx =t = *&(x) = —iS^(x). We consider the following special 
form of the WT relation (25): 



iti(Q\55C tot (x)\Q) = J dV (Q\T[5C tot (x')5C tot (x)]\Q) . (28) 
Here, 5C tot (x) and 55£ tot (x) are defined as 

5C tot (x) = 5C £ (x) = i[N(t), C £ (x)] tx=t = i(ee)v [*(x) - &(x)] (29) 

and 

55C tot = 55C £ = i[N(t), SC £ (x)} tx=t = -(ee)v [v(x) + &(x)] . (30) 
Thus, the WT relation is written in the form of 
— (ee)v 



v = 



2h 

(ee)v 



J dV [G n (x - x 1 ) + G 22 (x - x') - G 12 (x - x') - G 21 (x - x')\ 
[G 11 (p = 0)-G 12 (p = 0)] , (31) 



h 

where we have used the relations (16) in the last line 
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2.3 Diagonalization of unperturbed Hamiltonian and the Ward-Takahashi re- 
lation at tree level 



In this subsection, we investigate the system at tree level. First let us determine 
the order parameter at tree level. Form the unperturbed Hamiltonian (12) and 
the first term in the perturbative Hamiltonian (13), we have the diagram at 
tree (zero-loop) level for the condition (6), corresponding to the expression 

' (n\(p(x)\n) 

where /x is the chemical potential at tree level. As this equation must be true 
for any x, /x is determined as 

Ho = gv 2 -ee. (33) 



i J dV G(x - x') (-fio + gv 2 - ee) = . 



(32) 



Next, we construct the quasi-particle picture at tree level by diagonalizing 
the unperturbed Hamiltonian (12). The unperturbed quantum field ip(x) is 
expanded in terms of annihilation-operators a p (t) as 



2 °° 

<f(x) = 7p E a p {t)e^ p 



(34) 



Here, V is finite and the periodic boundary conditions on field operators (2) 
restrict p as 



P 



2nh 



\Tlxi fly, fT'Z ) 1 



n x ,n y ,n z : integers . 



(35) 



The symbol J^p^-oo stands for X^., ny ,n z =-oo which includes n x = n y = n z = 0. 
We emphasize that if one employs the Bogoliubov's prescription, the quantum 
field does not contain a (t) and breaks the CCRs, but our expansion contains 
it and the CCRs are held. The operators a p (t) and a p (t) are subject to the 
commutation relations of creation- and annihilation-operators, which follow 
from the CCRs for quantum fields (9) and (10): [a p (t) , aL, (t)} = <5 PP ' and 
other commutations vanish. The temporal development of tp(x) and a p (t) is 
generated by the unperturbed Hamiltonian (12). Substituting the expansion 
(34) and its Hermite conjugate into the unperturbed Hamiltonian (12), we 
obtain 



(e p + gv 2 + ee) a p a p + (a p a_ p + a p al p ) 



(36) 
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where we use the relation (33) and the notation of 
As is well known, the Bogoliubov transformation 



a p = u p b p + v p bi p , (38) 
where 



ILpi Vp i 

and 



e p + gv 2 + ee 1 
\ 2E P ± 2 (39) 



E p = yf* + 2gv 2 e p + 2 (e p + gv 2 ) ee + (eef , (40) 
diagonalizes the unperturbative Hamiltonian as 

oo 

H = E p b\,b p + const. (41) 

p=— oo 

We mention again that a crucial point here is that the mode with p = is 
included in the transformation. The infinitesimal symmetry breaking term (7) 
makes the transformation between a and b well-defined as long as e is kept 
finite. 

The field operator (p(x) is rewritten in terms of the creation- and annihilation- 
operators bp and b p as 



<p{x) = -= £ [b p u p er>- x ' E ^ +b p v p e' L n(p-^ p t) 

V y p=— oo 



(42) 



One can easily check that this quantum field and its Hermite conjugate satisfy 
the OCRs (9) and (10). 

The unperturbed matrix propagator of the field tp(x) on the vacuum |f2o), 
which is defined by the relation b p \Qo) = 0, is given by 



G (x-x') 



G ,u(x - x') G ,i 2 (x - x') N 

G Qt 2l(x-x') G^22{X-X') J 
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/ 



\ 



-i(n \T[ (p (x)^(x')]\n ) -i(Q \T[ v (xMx'm ) 

K -i(Q \T[^(x)^(x')]\Q ) -i(n \T[^(x)<p(x')]\n ) ) 
Its Fourier transform is defined by 



.(43) 



Go( P ) 



G ,u(p) G 0jl2 (p) 



d x 
{2irh)- 



-G Q {x)e^ x -^\ 



(44) 



and the explicit forms of its matrix elements are 



G f o,n(p) = Go,22(-p) = 
G f o,i2(p) = G r o,2i(p) = 



ut 



uj — uj p + \5 uj + ujp — i5 ' 



uj — Up + 18 u + Up — iS ' 
where ^u; p = E p and <5 is an infinitesimal positive parameter. 



(45) 
(46) 



Next, let us check the WT relation (31) at tree level, using the propagators 
(45) and (46). The right-hand side (RHS) of the relation (31) is manipulated 

as 



(ee v 



[Go,ii(p = 0)-G 0) i2(p = 0)] 



(ee v 



1 



gv 



ee 



h \Up = Q J y Tiuj p= o 
2gv 3 (ee) + (eefv 



h 2 u 2 p =0 



gv 



hu. 



p=0 / 



(47) 



where the last equality comes from the quasi-particle energy (40). This way 
the WT relation (31) at tree level is confirmed. We observe that the contri- 
bution of the zero-energy mode (the NG mode) is vital. In other words, for 
a finite volume system with a discrete spectrum, the Bogoliubov's prescrip- 
tion without operators representing the NG mode can not preserve the WT 
relations. 



2.4 Unitarily inequivalent vacua 



In the previous subsection, we considered the vacuum \Q ) associated with the 
operator b p as a physical one at tree level. There is another vacuum, denoted 
by |0), which is annihilated by a p : a p \0) = 0. We evaluate the inner product 
(0|f2 ) below, in the limit e — > but keeping finite V. 
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For this purpose, we will give explicit transformation to relate the two vacua, 
\Q ) an d |0), in Appendix A. We recapitulate Eq. (A. 43) which is the conclu- 
tion of Appendix A: 



|^o) = 



cxp 



■y OO ' 
p= — OO 



cxp 



1 °° V 

1 x - Up 4. 

2 £ 

z p=-oo U P 



p t t 



|0>, 



(48 



where the symbol Z)^L-oo' m eans summation without p = 0. Under the limit 
£ — > 0, u p and t> p with p ^ are finite, vq/uq becomes 1 and uq is divergent 
as s~i. Then one finds that 







0) 



(49) 



A conclusion in this subsection is that the vacua an d |0) are orthogonal 
to each other in the limit e — > even for finite volume system. This means 
that the Fock space built on |0) is unitarily inequivalent to one built on |f2o): a 
state ^ n |^ ) 

is never obtained by the superposition of a p ™'|0), where n and ml 
are integers. We emphasize that the existence of inequivalent representations 
does not necessarily require the thermodynamic (infinite volume) limit in this 
model. This clearly shows that the choice of the unperturbed Hamiltonian and 
the associated vacuum is essential for the determination of physical quantities 
even in a finite Bose-Einstein condensed systems. For example, one can get the 
Bogoliubov spectrum only for the vacuum in a broken phase |fi ), but never 
for that in a symmetric phase |0). The fact that the Bogoliubov spectrum and 
the Bogoliubov transformation are observed [13,14] forces us to choose the 
Bogoliubov vacuum |f2 ). 



2.5 The Hug enholtz- Pines theorem at zero temperature 



In this subsection, we derive the HP theorem from the WT relation at each 
loop level. First, we define Sfj, as the quantum correction to /xo (the chemical 
potential at tree level), 



ji — jiQ -\- 5 ji — gv 2 — ee + 5jx , 



(50) 



and Sfi will be determined later. We define the unperturbed Hamiltonian with 



H = J d 3 x tp\x) (-|^V 2 - frj <P(x) 
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2 

H — — \Alp\x)lp(x) + ip(x)ip(x) + (p\x)<p\x)} 



while the perturbative Hamiltonian H mt includes the 5[i term, 

H mt = v (-fi -5/j, + gv 2 - eej Jd 3 x + y? f (:r)} 

+ 9 V J d 3 ^ </7^(x)</7^(x)</?(x) + </?^(x)</7(x)</?(a;)j 
+ ^ J d 3 x (p^ (x)(p^ (x)(p(x)(p(x) 



(51) 



y d 3 x5fnp J< (x)(p(x) . 



(52) 



Since the unperturbed Hamiltonian H is not changed from that at tree level, 
one can adopt the quasi-particle picture and the unperturbed propagator at 
tree level, given in Subsection 2.3. 

Consider the Schwinger-Dyson equation, 



G-\p) = G 1 (p)-Z(p) + 5 -£l, 



(53) 



where G(p) and are the full propagator and the self-energy from the loop 
diagrams, respectively, and / is a unity 2x2 matrix. Both of G(p) and are 
matrices and can be calculated with the elements of the unperturbed matrix 
propagator (45) and (46) and the perturbative Hamiltonian (52). The matrix 
elements of have the properties of 



En(p) = E 22 (-p), Ei 2 (p) = S 2 i(p) , 
where the matrix form of the self-energy is written as 



(54) 



E(p) = 



£n(p) 2i 2 (p) 



^S 2 i(p) S 22 (p) 



(55) 



Let us rewrite the RHS of the WT relation (31) at any loop level, we obtain 
that 



{ee)v 



[G n (p = 0)-G 12 (p = 0)} 



(ee)v G^ u (p = 0) + G j 2 (p = 0) - E n (p = 0) - E 12 (p = 0) + f 



ft 



|G" 1 (p = 0)| 



(56) 
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where \G 1 (p)\ represents the determinant of G 1 (p)- Since Go has the form 
in Eqs. (45) and (46), one can obtain the explicit form of Gq 1 (p) as 



G \P) 
( 



(uO - UJ p )Up — (u + UJ p )Vp (u + UJ p )UpV p — (u — UJp)UpVp 



y(0J + UJp)u p Vp — (u — UOp)UpVp (uo - UO p )v 2 — (u + Up) 



• (57) 



u; 



The determinant of G (p) is calculated as 



G-\P) 



-uo — uo. 



( 2 2 \ - , / s. S JJj 

UO - UOp{Up + Vp) - Sn(p) + — 

- [2uo p UpVp - £i 2 (p)f 

-uo 2 + [S u (p) - E 22 (p)] + e 2 + 2^ 2 e p 



■ p ( u J + w J)-E 22 ( P ) + ^ 



+ -f + 



gv 



h h 



25 n 



+ En(p) + S 22 (p) 



2^ 
ft 



£12 (p) 



+■ 



ft 



2#v 2 25/i 



ft 



ft 



+ Sn(p) + S 22 (p) 



+ 



ft 2 



(58) 



We dropped the quadratic terms of the self-energy and Sfi in the second line 
of Eq. (58), because they are quantities of higher order. This way the RHS of 
the WT relation (31) at any loop level is organized as 



(ee)v 



[G n (p = 0)-G 12 (p = 0)] 



(ee)v j 2gv 2 5/j, 
ft ~ ~h 

2gv 2 \ 6 (A 



+ E n (p = 0) + E 12 (p = 0) 



(ee) 2 v 



ft 2 



x 



+2 



r y -f + £n(p = o)-£ 12 (p = o) 

£e f gv 2 5a ^ . .1 (ee) 2 



(59) 



The WT relation (31) requires that this should be equal to v, or equivalently 



(2gv 2 + ee) (-^ + S n (p = 0) - S 12 (p = 0)^ = 



(60) 



Since the nonvanishing order parameter v ^ in the limit e — * is under 
consideration, we finally attain the HP theorem 
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-P = E n (p = 0) - E 12 (p = 0) 



(61) 



We again stress that the HP theorem is proven from the WT relation in the 
e — > limit with finite V here, not in the thermodynamic (infinite volume) 
limit. 



2.6 The Hug enholtz- Pines theorem at one-Loop level 

In this subsection, let us check the HP theorem at one-loop level explicitly. 
According to the previous subsection, the HP theorem is a result of the WT 
relation. The loop expansion keeps the WT relation at each loop level. So, we 
naturally expect the HP theorem at one-loop level. 

From the condition (6) at one-loop level, S/i is given by 

dn = 2igG 0:11 (x-x)+igG ,i 2 (x-x) = — — . (62) 

V p=-oo ltj P 

The elements of the matrix self-energy at one-loop level are obtained as 



S n (p) = E 22 (-p) 



^ y p'=-oo 2cUp/CJ p _p/ 



/(icy, a; p _p/) 
- cjp' — ^p-p' + i5 



+ 



E 



S i2 (p) = S 2 i(p) 



2„,2 oo 



E 



h V p / = _ 00 2cjp/cjp_p/ 
5 g ^ 



f(-U>p>, -Up-pf) 

(63) 

/l (ay , CUp_p' ) ft- ( — iOpf , — CUp-p/ ) 

a; — u p i — ujp-p' + i5 cu + ky + ^p-p' — i<5 

(64) 



TiV „, m 2uy ' 

p'=— oo P 

where we have introduced functions f(uj p /,ujp_pi) and ft(uy', uy_ p ') as 



f (iVp> , Up-p>) 



3epiep.pi gv- 

72 Up'Up-p' -+- -T~2~{£p' + e p-p') n ^2 

/ i 



g 2 v 4 



Tr 



9 V i \ 1 / ■> 
^v^V + ^p-p') + j\ e p lUJ p-p' + e p-p'^p') 
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+ 



st 1 3e p ' 36p_p/ 



ft (ay, ujp-pi) — 



h \ h 

ft 2 
ft V ft 



+ 



ft 



+ 



2^ 
ft 



+ UV + id 



p-P 



+ 



3 fee) 



h 



2 ' 



zcUp'O'p-p' H — -y- 



ft^ 



-p—p' 



-\2 



+ 



2(ee) 
ft 2 



(65) 



(66) 



Then, we have: 



£n(p = 0) 

«?V ~ / 1 \ (6e 2 , rt „ 4^ 2 e p , 2(?V ee /3ew 6(ee) 2 l 
S i2 (p = 0) 

00 gv 2 



^JL_ V f-. (68) 
The following quantity, appearing in the RHS of Eq. (61), is simplified as 



S n (p = 0) - E 12 (p = 0) 

g 2 v* ~ / 1 W2e 2 , 4566 P , 4£e^ 2(ee) 2 ^ 

g 2e p , + 3^ 2 + 2ee 

^ ^ _g_ ^ 26 P , + 3gv 2 + 2e6 



ft V p,t"oo <V ft V jyt^oo 2uy 
_. _0_ 2e p / + gv 2 + 2ee 

~ hV f ^ 2E V , 
p — — oo v 

= ~h' 

Thus, the HP theorem has been confirmed at one-loop level. 

We comment on a possible choice of the unperturbed Hamiltonian H' Q which is 
obtained from replacing /i in Eq. (51) with /i given in Eq. (50). Then instead 
of H int we have the perturbative Hamiltonian H' mt without the Sfj, term. One 



(69) 
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can repeat the same procedure of diagonalizing the unperturbed Hamiltonian 
as in Subsection 2.3. As a result, the unperturbed Hamiltonian is given by 







£ E' p b'X + const., 



(70) 



where 




\J (e P + gv 2 + 5fif - g 2 v 4 . 



(71) 



We find a difficulty because E' becomes complex in soft momentum region at 
one-loop level. Recalling e p = p 2 /2m and taking the limit of e — > 0, we have 
at p = 0, 



As is seen from Eq. (69), 5fi is determined to be positive for positive g, which 
implies that E' p=0 is pure imaginary. Thus the choice of the unperturbed 
Hamiltonian HL does not offer a consistent treatment. 



3 Finite temperature case 

We extend the discussions on the WT relations and the HP theorem at zero 
temperature to finite temperature case. We employ the TFD formalism to 
describe equilibrium situations. As will be seen, the TFD formalism is suitable 
for our purpose, because it is formulated as a canonical operator formalism of 
quantum field. 

3.1 Thermofield dynamics 

In this subsection, we give a brief review of TFD formalism. In TFD, thermal 
degrees of freedom are introduced by doubling each degree of freedom through 
the tilde conjugation. Thus, with an operator A we associate its tilde conjugate 
A according to the tilde conjugation rules [4,18]: 




(72) 



(AB)~ 



AB, 



(73) 
(74) 
(75) 




(At)~ 
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(A)~ = aA, (76) 
\Q)p" = \Q)p, (77) 
p(n\~ = (%> (78) 

where c\ and C2 are complex c- numbers, and ^(Q] are the thermal vacua, 
and a = 1 for bosonic A, while er = —1 for fermionic A. The total Hamiltonian 
H is given in terms of the non-tilde and tilde Hamiltonians as 

H — H — H , (79) 



and the total Lagrangian L is also given as the non-tilde Lagrangian minus 
the tilde one, 

L — L — L. (80) 



The thermal average is expressed as the thermal vacuum expectation value in 
TFD: ^(filAlJl)^. For simplicity, consider one bosonic degree of freedom whose 
creation and annihilation operators are represented by a and a', respectively. 
When the density matrix p is expressed in the form of p = f a ^ a , we have 

^1^ = ^ (81) 

and 

n = p (n\a^a\n)p = J-j, (82) 



using the properties of thermal vacua. If / = exp(— f3E) where (3 is the inverse 
temperature and E is an energy, ^(QI^IQ) ^ represents the thermal expectation 
value at finite temperature. One can easily generalize this argument to multi- 
mode cases. 



For convenience for following calculation, let us introduce the matrix notation 
0f and 4>i which represent tilde and non-tilde quantum field f and dp as 



3, 



# = 4# = -<fc, (83) 
where 
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01 = V?, 02 = $ , 

4>\ = <p\fo = <p. (84) 

Then, the thermal propagators with four indices and with sixteen components 
are defined by 



G^ix-af) = -ip(n\T[^{x)^{x'm)^ (85) 

where /j,, v — 1, 2 and i,j = 1, 2. The Fourier transformed thermal propagators 
are also defined as 



G MW=[j0pGfa(x)ei<r^K (86) 



One can derive the following relations of the Fourier transformed thermal 
propagators by general properties of TFD: 

G%M = G% 22 {- P ) , G%M = Gfaip) , Gf(p = 0)=0. (87) 

3.2 The Ward-Takahashi relation at finite temperature 

In TFD, the total Lagrangian £ density is given by 

£ = £-£. (88) 

We consider following continuous infinitesimal transformations: 

ip(x) — > ip\x) = ip(x) + £5ip(x) 

$(x) -> tp'(x) = $(x) + t5$(x) . (89) 

The Nother theorem at finite temperature is obtained from the infinitesimal 
transformations (89) as 



d - 

—N°( y x) + V-N( y x) = S£, (90) 
at 



where 



d£ BC 

N»(x) = ————5^(x) - -^4 ^(x) (91) 

d(d^(x)) d(d^(x)) 
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and 



5C = 5C-5C. (92) 
Here, 5C was defined in Eq. (18) and SC is denned as 

£5£ = £[*'(x)] - £[*(*)] . (93) 



One can easily obtain the WT relation at finite temperature by replacing 5C 
with 5 1 in Eq. (25): 



J2 ihp{Q\T[^(x 1 ) ■ ■ ■ 5V(x a ) ■ ■ ■ V(x n )]\Q)p 

a=l 

J d 4 xp(Q\T[5C(x)^(x 1 ) ■ ■ ■ *(x n )]\n) p . (94) 



Now we consider the following WT relation with respect to the global phase 
transformation, i.e., S^(x) = i^(x) and 5^f(x) = —i^f(x) in Eq. (89). To deal 
with the spontaneous breakdown of the symmetry at finite temperature, we 
also introduce an infinitesimal symmetry breaking term with respect to the 
tilde field as 

C £ = (ee)v [^(x) + &(x)] , (95) 



which is a tilde conjugate of C £ (7). Our starting total Lagrangian density at 
finite temperature is given by 

Aot — Aot — Aot , (96) 



where £ to t was defined in Eq. (8) and C tot is defined as 

Act = I + 4 . (97) 



The Nother charge of the total Lagrangian £ tot (96) with respect to the global 
phase transformations is obtained: 



N(t) = N(t) - N(t) = - Jd 3 x {^(x)^(x) + ¥(x)V(x)} . (98) 
This generates the infinitesimal transformation, e.g., 
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5C tot = i[N(t),Ctot(x)] = i(ee)v [*(x) - &(x)] , 
55£ tot = i[N(t),5jC tot (x)\ = -(ee)v \*(x) + &(xj 



(99) 
(100) 



and 



5£ tot {x) = i[N(t),£ e (x)] = i(ee)v - &(x) + *(x) - ^(rr)] . (101) 



One can derive the following restricted relation from the WT relation at finite 
temperature (94), and from Eqs. (99), (100) and (101): 



thfs(n\55C tot (x)\n} p = j d 4 x' p (n\T[5jC tot (x')5jC tot (x)]\n) p . (102) 

We rewrite this relation in terms of the thermal propagators: 



(ee)v 



v = 



2h 



d 4 x' (79*11 (x — x') + G^22 ( x ~ x> ) ~~ G l ^i2(x — x') — G^2i(x 

+ G^nix — x ) + G^22( x — x ) — G^ 12 (x — x ) — G^2i(x — x ) 
(ee)v 



Gf, n (p = 0)-Gf A 2(p = 0) 



3. 3 The Ward- Takahashi relation at tree level in finite temperature case 

Let us introduce the operators {£ p , £„} and {£ p , £ P } related to the operators 
{b p , bp} and {b p , b p } by the following thermal Bogoliubov transformation: 



bp c p £p + s p^p 1 



(104) 



with 



e 2 



A 



(105) 



where the energy £" p is given in Eq. (40). Since c p — s p = 1, we find that 

[£pj w] = [Ipi w] = ^pp' an d the other commutation relations vanish. With 
these operators, the unperturbed Hamiltonian at the finite temperature in 
TFD is given as 



Ho= it E p (b% - b%) = £ E p (# p _ 

p= oo p=— oo 



(106) 
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The field operator (p(x) is rewritten in terms of the operators {£ p , £ P } and 
{i P , IJ} as 



1 _ 



+ (cp4 + splp) V~* <P ""~ £W ] • (107) 

We construct the unperturbed propagators with this field operator and its 
Hermite and tilde conjugate: 



= -i^nolT^x)^)]!^, (108) 
and its Fourier transformed one 



G 



s.«(p) = / ^^ ,(,) e i(— ) . (109) 



Here, the unperturbed thermal vacuum \ £lo)* is specified by £ p \Qo) ^ = £ p |f2 ) ^ = 
0. The explicit forms of the matrix elements, which are necessary to check the 
WT relation, are given as 



Gfc,ii(p)=G^ >22 (-p) 



c 2 s 2 \ 

P P \ 2 

1 u„ 



\u — Ld p + iS uj — ujp — iSJ p 
+ ( ( 11Q ) 

\U + LOp + 10 U + LOp — 10 J y 



c 2 s 2 
s> _l *p_ 



^ u — Up + 18 u + Up + iS 



u + uj p — iS u — Up — iS 
At p = they are reduced to 



2 2 \ 

CP SP ^UpVp. (Ill) 



Gf, , 11 ( P = 0)=Gl\ 22 (p = 0) 



1 



U!p=Q 

2 



(c p=0 s P =o) ( M p=0 + v p=o) 



u 



p=0 + v p=0 



u. 



p=0 



G , o,ii(P = 0) = Go,22b = 0), 



(112) 



23 



Gf J0 , 12 (p = 0)=G 1 ^ n (p = 0) 
2 



UJp=o 

2u p=0 v p=0 



2 2 
c p=0 — s p =o ) Up =0 V p= o 



= G 0l i2(p = 0) = Go,2i(p = 0). 



(113) 



Substituting Eqs. (112) and (113) into Eq. (103), we find that the WT relation 
at finite temperature holds at tree level: 



[ee)v 



GSn(p = 0)-Gg 12 (p = 0) 



[ee)v 

h 



[G n (p = 0) - G 12 (p = 0)] = v(lU) 



In other words, at tree level, the Goldstone theorem also holds at finite tem- 
perature. 



3.4 The Hug enholtz- Pines theorem at finite temperature 



In Subsection 2.5, the HP theorem at T = was derived from the WT relation. 
Now we will generalize this story to finite temperature case. We considered 
the quantum corrections of the chemical potential (50) at T = 0. At finite 
temperature one has to take account of thermal corrections as well as quan- 
tum ones. Explicitly Sfj, in Eq. (50) now represents the thermal and quantum 
corrections to the chemical potential at tree level. 

We consider the following Schwinger-Dyson equation at finite temperature: 



G,\P) = G^ip) - ^{p) + ^J 4X4 , (115) 

where Gpip) and S,g(p), having the indices (/x, v) and (?, j) are the full prop- 
agator and self-energy at finite temperature, respectively, and J 4X 4 is a unity 
4x4 matrix. In TFD, one can find the following relations in general, 



Gf(p = 0) = Gf(p = 0) , Gf(p = 0) = Gf(p = 0) = , 
Ef{p = 0) = Ef(p = 0) , 4 2 (p = 0) = Hf{p = 0) = , 
Gj} u (p = 0) = Gf 22 (p = 0) , Gl\ 2 (p = 0) = Gf 21 (p = 0) , 
EjJ u (p = 0) = Z£ j22 (p = 0) , Eg 12 (p = 0) = Zf i21 (p = 0) . (116) 



Let us calculate the RHS of Eq. (103), 
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h 




(ee)v 




h 




(ee)v 




h 



{<V> = o)} - {<v> = °)} J W> = °) + O = o) 

{g,Z(p = 0) - Eg u (p = 0) + - {G^ > = 0) - £g 12 (p = 0)f 



x 



G-,Z(P = 0) + = 0) - E£ u (p = 0) - Sg 12 (p = 0) + ^ 



From Eqs. (112) and (113), we find that this expression is very similar in form 
to that at T = in Eq. (56), except that the self-energy terms Sy u (p = 0) 
and S^ 12 (p = 0) replace £n(p = 0) and S 12 (p = 0) there. Thus the WT 
relation at finite temperature is equivalent to the following HP theorem at 
finite temperature: 



(117) 



5/i 



E^n(p = 0)-E^ 12 (p = 0), 



which is very similarly as Eq. (61) has been derived at zero-temperature. 



(118) 



3.5 The Hug enholtz- Pines theorem at one-loop level infinite temperature case 



Finally, in an explicit calculation at one-loop level in finite temperature case, 
we check the conclusion in the previous subsection, i.e., the HP theorem or the 
WT relation. The condition (6) fixes Sfj, at one-loop level in finite temperature 
case (See Eq. (62) at T = 0), 



5fi = 2igG/3 A11 (x - x) + ^^,0,12(2; - x) 

gv 2 + 2ee 




(119) 



To investigate the HP theorem, it is necessary to have only the two matrix 

-111 

J /3,l 



elements of the self-energy, Si 1 11 (p) and Y$ 12 (p). They read as 



a 2 ?; 2 00 



ft Z p^oo V 2^p'CJ p _p/ 



/ 2 2 2 2 \ 

I C p C p— p' S p S p—p J 

, f 2 2 22 \ I 

' \ c p s p—p s p ,c p _ p ,j < 



p-p 



f(-up 



-00. 



P-P 



00 — Up 



t^p-p + iS 00 + cup + tu p -p — i5 
f (up , 0J p _p ) /(— uy , —ojp-p ) 1 



00 — oup + ojp^p + i5 



U) + 0i) z 



ojp_p — i5 
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6p ' + ^ 2 + £e (120) 



n 2 2 oo i 



h((jjp> , UJp-pi) h(—u p >,—LJ p _ p > 



fl pl= _ OD V liOpiLOp-pi 
( 2 c 2 _ g 2 g 2 \ J 

V p p-p p p-p / | ^ _ ^ _ ujp-j/ + i5 uj + ay + cjp_ p / — i<5 J 

/2 g 2 _ s 2 c 2 )[ h(uj p ,,uj p _ pl ) h(-LO p/ , -UJ p - p ') 

V p p-p p p-p / I _ _|_ _|_ j(5 ^ _|_ _ uj p _ pl — iS 
(3E P ,\ gv 2 




Using c p — s p = 1 and referring to the calculations in Subsection 2.6, we can 
easily find that the RHS of Eq. (118) is given by 



h.^V 2 2E P / 

p— — oo \ / P 



(122) 



Thus it has been confirmed that the HP theorem or the WT relation is satisfied 
at one-loop level. In other words, our approximate calculation scheme respects 
the Goldstone theorem at the one-loop level in finite temperature case. 



4 Summary 



We constructed an approximate scheme satisfying the Goldstone theorem for 
Bose-Einstein condensed gas in a box with the periodic boundary condition 
and generalized it to finite temperature case using TFD. When energy spec- 
trum is continuous as in homogeneous situation, one takes the thermodynamic 
limit and then the Bogoliubov's prescription do not cause serious problem for 
the Goldstone theorem, since the NG mode, which is p = mode, is just a 
point in momentum integral and omitting it does not affect the result. On the 
other hands, when a system has a finite volume so that energy spectrum is 
discrete and that the concept of the thermodynamic limit is not applicable, we 
found that it is necessary to deal with the NG mode explicitly for satisfying 
the WT relation, that is, the Goldstone theorem. Without the thermodynamic 
limit, the HP theorem is also derived from the WT relation. Furthermore the 
NG mode contributes essentially to the unitarily inequivalence of vacua in the 



(121) 
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case of the SSB for finite volume system. Then it is found that the scheme for 
the trapped dilute Bose system satisfying the Goldstone theorem also need 
explicit treatment of the NG mode. 
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A Bogoliubov transformation and vacua 

In this appendix, we give explicit relation between |0) and \Q ) which are the 
vacuum of a p and b p respectively [19]. 

The Bogoliubov transformation (38) is written as the following matrix form: 




(A.l) 



Next, we introduce an operator U p which induces the transformation (A.l): 
b p = U p a p Ul , bl p = UpalpUl , (A.2) 

where U p must be unitary, 

UpUl = U P U P = 1 . (A.3) 
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Here, 1 is the identity operator. Eqs. (A. 2) imply 
b p U p \0) = 0, <0|C/+6l„ = 0. 



(A.4) 



We define the unitary operator U as 

oo 

U= ]JU P , 

p=0 



(A.5) 



where Il^Lo contains the suffix p = 0, and U is obviously unitary. Then the 
relation between \Q ) and |0) is given as 



\n ) = u\o). 



(A.6) 



First, we consider the p ^ part of U. To determine an explicit representation 
of Up in terms of a p and a^ p , we define the coherent state of a p and a_ p , 



\z(p)) = exp[z p a P ] exp[z_ p al p ]|0) 



(yZ — p(l_p 



|0),(A.7) 



where z p and ^_ p are eigenvalues of a p and a_ p , i.e., 
a P \ z (p)} = z p\ z (p)) , a - P \ z (p)} = Z - P \ Z (P)) ■ 



(A.8) 



The coherent state has the properties of 

a U z (p)) = 9 Zp \z(p)) , oL p \z(p)) = d z _ p \z(p)) , 
(z(p)\a P = (z(p)\z* p , (z(p)\al p = (z(p)\z*_ p , 
{z(p)\a p = d z *(z(p)\ , (z(p)\a- p = d z *_ (z(p)\ , 



and 



<0|z(p)) = l. 



(A.9) 
(A.10) 
(A.H) 



(A.12) 



The completeness condition of the coherent state reads 

Jdfi(z(p))\z(p))(z(p)\ = j&n{z{p)) l\z p )(zp\®\z_p)(z_p\] = l, (A.13) 



where the measure is defined as 
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d/i(z(p)) 



-\Z V \ 2 dX P d V P 



7T 



e -\*-P 



2 dx_ p dy. 



IT 



z± P = x± p + iy± p . 



(A.14) 
(A.15) 



We also mention the following formula of a generalized Gaussian integral: 



J dn(z(p)) exp -- (z T {p)Mz{p) + z* T {p)N*z*{p)) + u T z{p) 

: [det(J - MiV*)P 
x exp 

+ v* T (I - MN*)- 1 !! 



v* L z*(p) 



■- u T (I - N*M)~ 1 N*u - ^v* T (I - MX') '.!//•• 



(A.16) 



Here, M and N* are 2x2 symmetric matrices, 



M T = M , N* T = N* , 



(A.17) 



7 is a unity 2x2 matrix, and z, u and v are c-number 2-components vectors, 

\ 



z(p) 



u 



Z-, 



u_ 



(A.18) 



Now, we are led to the operator U p explicitly in terms of a p and al p . One can 
easily derive the following relations from Eqs. (A.2)-(A.4) and (A. 8): 



(b P - z p )U p \z{P)) = , (bl P ~ d z _ p )U p \z{p)) = . 



(A.19) 



Putting the expression (A.l) into Eqs. (A.19), we obtain 

(u p a p - v p aL p - z p )U p \z(p)) — , 
(-v p a p + WpaL p - d z _ p )U p \z(p)) = , 

whose inner product with the coherent state (w(p)\ are given by 



(A.20) 
(A.21) 



(u p d w * - v p w*_ p - z p )(w(p)\Up\z(p)} =0 , 
{-v p d w * + u p w*_ p - d z _ p )(w(p)\U p \z(p)) =0 . 

One finds that the following solution is admitted, 



(A.22) 
(A.23) 
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(w(p)\U p \z(p)) 

— ex P \lw* T (p)X p w*(p)-lz T (p)X p z(p)+w* T (p)Y p z(p)] , (A.24) 



where 





(A.25) 



and 



w(p) 



( 



W-p 



w. 



V 



) 



(A.26) 



In the definition of the matrices X p in Eq. (A.25), we used the Pauli matrix: 



The w- and ^-independent normalization factor in Eq. (A.24) is fixed by the 
relation, 



where the formula (A. 16) can be applied to the last expression of the integral. 
Thus we obtain an explicit operator representation of U p with the help of the 
relation A(w* p , w*_ p , z p , z_ p ) = (w(p)\:A(a\ t ,al. p ,a p ,a- p ):\z(p))/(w(p)\z(p)) 
where the symbol : • • • : implies the normal ordering to keep the creation 
operators to the left of the annihilation operators, 




(A.27) 



exp (w*(p)w'(p)) = (w(p)\w'(p)) 

= {w(p)\U p Ul\w'(p)) 



J Mz{p)) (w(p)\U p \z(p))(z(p)\Ul\w'(p)) 

J Mz{p)) (w(p)\U p \z(p)) ((w'(p)\U p \z(p))Y ,(A.28) 



U p = — : exp ia fl \p)X p a ] \p) 



±a T (p)X p a(p) + a^(p)Y p a(p)j : , (A.29) 



where 




(A.30) 
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Next, we consider the p = part of U. The coherent state of the p = mode 

as 



\z(0)) = exp[z o al]\0) 



E 

1=0 



(zoaiy 



|0>, 



(A.31) 



where 



and 



a |z(0)) =Zq\z{Q)) , 4\z(0)) = d Zp \z(0)) , 
(z(0)\a o = d z *(z(0)\, (z(0)\4 = z* o (z(0)\ , 

\z(0)) = 1. 



(A.32) 
(A.33) 

(A.34) 



The completeness condition of the coherent state reads 

Jd»(z(o)) = i, 



(A.35) 



where the measure is defined as 



d^(0))=e-l-l 2 ^^, 

7T 

z = x + iy ■ 
The following Gaussian integral is relevant here: 



(A.36) 
(A.37) 



J dn(z(0)) exp \~ (M'zl + N'*z* 2 ) + u'z + v'*z* 



= (1 - M' AH exp 



1 N'*u' 2 + M'v 



*2 



2u'v' 



1 - M'N' 



(A.38) 



where N', M'*, u', v'* are constants. In a manner similar to the manipulations 
of the p 7^ part of the U, we obtain the following equations, 



(u d w * -v Wq -zo)(w(0)\U o \z(0))=0, 
(-v d w * + u w* - d zo ){w(0)\U o \z{0)) = . 

One easily find the matrix elements of U as 



(w(0)\U \z(0)) 



exp 



2 M «0 



(A.39) 
(A.40) 



(A.41) 
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so the representation of Uo in terms of a p and a is given as 



U = 



: exp 



1 ^0 / +2 2 \ , 1 t 

- — (a -a ) + —a a 

2u u 



(A.42) 



As a summary of this appendix, the vacuum of the quasi-particle, described 
by the b p operators including the p = case, is explicitly given by 

\too) = U\0) 




where l\ P < ' = _ OD ' and SplLo' mean product and summation without p = 0, re- 
spectively. 
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